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ABSTRACT 

We present the analysis of the spectroscopic and photometric catalogues of 11 X- 
ray luminous clusters at 0.07 < z < 0.16 from the Las Campanas / Anglo- Australian 
Telescope Rich Cluster Survey. Our spectroscopic dataset consists of over 1600 galaxy 
cluster members, of which two thirds are outside ^200- These spectra allow us to assign 
cluster membership using a detailed mass model and expand on our previous work on 
the cluster colour-magnitude relation where membership was inferred statistically. We 
confirm that the modal colours of galaxies on the colour magnitude relation become 
progressively bluer with increasing radius d{B — R)/drp — —0.011 ± 0.003 and with 
decreasing local galaxy density d{B — R) / dlog{Yj) = —0.062 ± 0.009. Interpreted as an 
age effect, we hypothesize that these trends in galaxy colour should be reflected in mean 
H(5 equivalent width. We confirm that passive galaxies in the cluster increase in H(5 line 
strength as dRS/drp = 0.35 ± 0.06. Therefore those galaxies in the cluster outskirts 
may have younger luminosity- weighted stellar populations; up to 3 Gyr younger than 
those in the cluster centre assuming d[B — R)/dt = 0.03 mag per Gyr (Kodama & 
Arimoto 1997). A variation of star formation rate, as measured by [Oii]A3727A, with 
increasing local density of the environment is discernible and is shown to be in broad 
agreement with previous studies from the 2dF Galaxy Redshift Survey and the Sloan 
Digital Sky Survey. We divide our spectra into a variety of types based upon the 
MORPHs classification scheme. We find that clusters at z '--^ 0.1 are less active than 
their higher redshift analogues: about 60 per cent of the cluster galaxy population is 
nonstarforming, with a further 20 per cent in the post-starburst class and 20 per cent 
in the currently active class, demonstrating that evolution is visible within the past 
2-3 Gyr. We also investigate unusual populations of blue and very red nonstarforming 
galaxies and we suggest that the former are likely to be the progenitors of galaxies 
which will lie on the colour-magnitude relation, while the colours of the latter possibly 
reflect dust reddening. We show that the cluster galaxies at large radii consist of both 
backsplash ones and those that are infalling to the cluster for the first time. We make 
a comparison to the field population at z ~ 0.1 and examine the broad differences 
between the two populations. Individually, the clusters show significant variation in 
their galaxy populations which we suggest reflects their recent infall histories. 

Key words: surveys - catalogues - galaxies: clusters: general - cosmology: observa- 
tions ~ galaxies: evolution 
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Clusters of galaxies represent an ideal cosmological labo- 
ratory. They allow the study of large numbers of galaxies 
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10"^), all at common distances. Since they are visible out 
to high redshifts, they are also excellent tools for studying 
galaxy evolution and the role that environment plays in this 
evolution. Indeed, the environment of a galaxy is likely to 
be one of the most important factors that affect its star for- 
mation rate (Osterbrock 1960; Dressier, Thompson & Sheet- 
man 1985; Abraham et al. 1996; Balogh et al. 1997; Balogh 
et al. 1998; Small et al. 1998; Poggianti et al. 1999; Kodama 
et al. 2001; Lewis et al. 2002; Pimbblet et al. 2002 [P02 
herein]; Gomez et al. 2003; Balogh et al. 2004). However, it 
could also be that initial conditions are just as important 
as environment in determining the star formation history of 
galaxies. 

It is well-established that galaxies that reside in the 
high (iousity cores of clusters generally exhibit lower cur- 
rent star formation rates than equivalent field galaxies. 
The morphology-density relationship (T — S; Dressier 1980; 
Whitmore, Gilmoro & Jones 1993; Dressier ct al. 1997; 
Smith et al. 2005; see also Oomlcr 1974) can explain some of 
this trend. Broadly, it shows that latc-type galaxies, which 
are expected to correspond to star-forming systems, prefer- 
entially occupy low-density regions whilst conversely, early- 
type passive galaxies tend to be situated in regions of high 
density such as at the centre of galaxy clusters. The fraction 
of early types is also found to be correlated with global clus- 
ter structure: more relaxed clusters have a higher early-type 
fraction (Oemler 1974). The origin of the T — E relationship, 
however, has been the subject of debate for many years and 
it is also known that the morphologies of galaxies and their 
current star formation rates do not necessarily follow the 
same correlations in different environments (Balogh et al. 
1997; Poggianti et al. 1999; Couch et al. 2001; Lewis et al. 
2002; Gomez et al. 2003; Trou et al. 2003; Hogg et al. 2004; 
Christlein & Zabludoff 2005; see also Girardi et al. 2003). 

One obstacle in our understanding of the role of envi- 
ronment in galaxy evolution is that many studios to date 
have probed the cluster regime in complete isolation from 
their infall regions and the field at largo. Whilst there arc 
a few notable exceptions (e.g. Fasano et al. 2005; Rines et 
al. 2005; Christlein & Zabludoff 2005; Gerken et al. 2004; 
Abraham ct al. 1996), the lack of data at large radii from 
the centre of clusters represents a very serious impediment 
in our understanding of the role that environment plays. 
This work analyzes one dataset that is capable of simulta- 
neously probing the very high density cores of galaxy clus- 
ters, their infall regions (which typically have a projected 
density corresponding to a poor group) and out into the 
field: The Las Campanas Observatory and Anglo-Australian 
Telescope Rich Cluster Survey (LARCS; e.g. Pimbblet et al. 
2001; POl herein). 

LARCS is a long-term project to study a statistically- 
reliable sample of 21 of the most luminous X-ray clusters at 
z = 0.07-0.16 in the southern hemisphere*. The clusters are 
a random subsample of the X-ray brightest Abell clusters 
(Abell 1958; Abell, Corwin & Olowin 1989) selected using a 
cut of Lx > 3.7 X lO''* erg/s from Ebeling et al.'s catalogue 




Radius / Mpc 



Figure 1. Fraction of projected number density of LARCS galax- 
ies down to M* + 3 which arc cluster members as a fraction of the 
total projected number count. At large radii, > 3 Mpc, we can 
expect that as few as one in ten galaxies observed will belong to 
the cluster. LARCS allocates a large fraction (up to 70 per cent 
per pointing) of 2dF fibres to galaxies in the cluster outskirts to 
ensure that we obtain a reliable sample of galaxies at these radii. 



(1996). X-ray luminosity, Lx, is a good proxy for mass, so we 
should be selecting only the most massive clusters at these 
redshifts. However, we note that our sample may be contam- 
inated with less massive clusters that possess boosted Lx as 
a result of dynamical activity or the merging of subclusters. 
Although this range corresponds to only quite modest look- 
back times, it spans a 1.3 — 2.6 Gyr^^ period (or similarly, a 
1.0 - 2.0 Gyr period for a ft = 0.7, A = 0.7 flat cosmology) 
that has remained relatively under-explored, giving rise to a 
"gap" in the tracking of rich cluster evolution back in time. 
Such a long baseline is necessary if we are to catch clusters 
in their important phases, such as mid-merger, and having 
a sample of > 10 clusters will guarantee that we fairly sam- 
ple most of the important phases in their evolution (e.g. 
mergers, Lacey & Cole 1993). Accordingly, we are mapping 
the photometric, spectroscopic and dynamical properties of 
galaxies in rich cluster environments at 2: ~ 0.1, tracing the 
variation in these properties from the high-density cluster 
cores out into the surrounding low-density field beyond the 
turn-around radius. For the most massive clusters at 2; ~ 0.1, 
the turn-around radius corresponds to roughly 1 degree or a 
10 Mpc radius (O'Hely ct al. 1998; P02) and therefore we 
have obtained panoramic CCD imaging covering 2-degree 
diameter fields, as well as spectroscopic coverage of these 
fields (e.g. POl; O'Hely 2000; Pimbblet 2001). The imag- 



* Their declinations are constrained to be < +10 so as to be ac- 
cessible from both Las Campanas Observatory and Siding Spring 
Observatory. A reddening limit oi Ag < 0.3 mag is also imposed 
(based on values from Burstein & Heiles 1984). 



t Throughout this work values of Hq = 50 km s~'^Mpc~^ and 
qo = 0.5 have been adopted. Further, all quoted coordinates are 
J2000 compliant. 
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ing comes from B and _R-band mosaics taken with the 1-ni 
Swope telescope at Las Campanas Observatory, while the 
spectroscopy comes from the subsequent follow-up with the 
400-fibre 2dF multi-object spectrograph on the 3.9-m Anglo- 
Australian Telescope (A AT). 

The LARCS dataset permits studies into the role of 
environment and redshift evolution of cluster galaxies as it 
possesses the requisite wide-field observations coupled with 
dense sampling of individual clusters to effectively bridge the 
gap between cluster centres and their low-density outskirts 
at large radii whilst simultaneously providing the lower red- 
shift analogues to high redshift cluster studies (e.g. Pog- 
gianti ct al. 1999). It is vital to have this wide-field coverage, 
as an important constraint on the nature of the processes 
that decrease the star formation rates of infalling galaxies 
is the star formation efficiency in low density environments, 
of ~few galaxies Mpc~^ (Gunn & Gott 1972; Larson, Tins- 
Icy & Caldwell 1980; Barnes & Hernquist 1991; Moore et 
al. 1996; Quilis, Moore & Bower 2000; Kodama ct al. 2001; 
Lewis et al. 2002; Balogh et al. 2002; Bekki, Couch & Shioya 
2002; Gomez et al. 2003; Balogh et al. 2004; Christlein & 
Zabludoff 2004; Burstein et al. 2005). This density is typical 
of those galaxy groups that are just starting to fall into the 
cluster - i.e. at large radii from the cluster centre. Therefore, 
not only do we need to have panoramic coverage with both 
imaging and spectroscopy to encompass all environments 
within a cluster to understand the processes, but we also 
need a good rmmbor of highly sampled clusters themselves 
in order to trace all likely evolutionary pathways. LARCS 
is ideally situated to shed light on these areas due to its 
impressive spatial extent of up to 10 Mpc radius (see Fig- 
ure 1) that allows us to directly compare galaxies at very 
low density to those at very high density. We arc also in 
an excellent position to inter-compare clusters within our 
homogeneously selected sample to investigate the degree of 
uniformity amongst them. Moreover, by having a homoge- 
neous sample, we can co-add the individual clusters to try 
to improve the statistics of rare classes of galaxies in the 
ensemble population. Our sample will therefore provide a 
« ~ 0.1 benchmark to compare to higher redshift cluster 
samples. 

The format of this paper is as follows: in §2 our ob- 
servational procedure and data reduction methods are de- 
scribed. We address our selection function and derive com- 
pleteness limits in §3. We perform a dynamical analysis in 
§4 to determine cluster membership via a mass model and 
classify our galaxies into a number of spectral types. In §5, 
we examine the cluster colour-magnitude relation, analyze 
the environmental dependence of line strengths and present 
a breakdown of each cluster on spectral type. Our findings 
are discussed and summarized in §6. Appendix A presents 
our spectroscopic catalogues while Appendix B catalogues 
other large-scale structures and galaxy cluster candidates 
serendipitously found along the line of sight. 

2 THE DATA: OBSERVATIONS AND 
REDUCTION 

2.1 Optical Imaging at LCO 

High quality broad-band B- and i?-band CCD images of our 
clusters have been secured at Las Campanas Observatory 



(LCO) using the 1-m Swope Telescope. More details of the 
observations, reduction and analysis of these data are given 
in POl; here we briefly summarize the pertinent points. 

A 2 deg-diameter field around the cluster is imaged in 
21 over-lapping pointings to produce a mosaic of a region 
out to ~ 10 Mpc radius at the cluster redshifts. The images 
are reduced using standard tasks within iraf and are cat- 
alogued using the SExtractor package of Bertin & Arnouts 
(1996). We adopt mag_best from SExtractor as the estimate 
of the total magnitudes. To determine colours for the galax- 
ies we perform aperture photometry within 4" apertures, or 
~ lOkpc at the typical cluster redshift, on seeing-matched 
tiles using piiot within iraf. Photometric zeropoints are 
computed using the frequent observations of standard stars 
from Landolt (1992) interspersed throughout the science ob- 
servations. On photometric nights, the variation in the re- 
sultant colour and extinction terms are all within Icr of each 
other and the final photometric accuracy is better than 0.03 
mags. 

Finally, wo apply corrections for galactic reddening 
based on Schlegel et al. (1998). The final internal magnitude 
errors across the full mosaics are typically 0.03 magnitudes 
and always less than 0.06 mags (POl). The catalogue is ~ 80 
percent complete at a depth of i? ~ 22.0 and B ~ 23.0. 

Star/galaxy separation for the catalogue uses the ro- 
bust criteria described in POl: galaxies aie selected by re- 
quiring FWHM > 2.0" and CLASS_STAR < 0.1 from the SEx- 
tractor star-galaxy classifier. The resulting stellar contam- 
ination based on these criteria is estimated at < 3 percent 
(see POl). 

2.2 2dF Spectroscopy and Target Selection 

The spectroscopy used in the LARCS project comes from the 
2dF spectrograph (Gray & Taylor 1990)"^ mounted upon the 
AAT. Our 2dF observations are summarized in Table 1 and 
we present our catalogues in Appendix A (see also Pimbblet, 
Edge & Couch 2005). Briefly, these observations were per- 
formed over the course of ^12 clear nights from 1998 to 2002. 
We elected to use the 600V grating as it provides the neces- 
sary resolution to measure key line indices and precise veloc- 
ity measurement (~120kms~^). The restframc wavelength 
range with this grating is approximately 3700 5600A for our 
clusters, providing coverage of the major line indices of in- 
terest ([Oii]A3727A, [Oiii]A5007A, H5, H/3, Mgb, FeA5270). 
The mean seeing for these observations was ~1.5". Our ex- 
posure times arc typically 9ks per cluster per pointing in 
dark conditions with two pointings per cluster. These times 
are chosen in order to obtain a signal to noise ratio of 10- 
15 per pixel for our faintest galaxies {R ~ 20.5). Here, one 
pointing is defined as a combined set of 2dF observations, 
totalling 400 fibres. Our goal was to obtain two such point- 
ings per cluster in order to obtain approximately 150 galaxy 
cluster members with high quality redshifts - a number re- 
quired if we are to be able to identify small sub-clumps and 
rare populations such as E-|-A galaxies. 

The galaxies used as targets for the 2dF observations 

t The 2dF spectrograph is described in detail on part of the AAT 
website (http://www.aao.gov.au/2df/index.html) and in Lewis et 
al. (2002). 
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are chosen from the optical photometric catalogue derived 
from SExtractor. As the accuracy required by the 2dF fibre 
positioner is < 0.3", only galaxies whose SExtractor FLAG 
parameter (see Bertin & Arnouts, 1996) is less than 4 are 
used^; a higher value of the FLAG parameter would be a 
concern for the accuracy of the astrometric solution for that 
galaxy. 

With the aim of efficiently sampling the whole clus- 
ter population, the galaxy sample is split into four parts 
based upon radial distance from the cluster core and magni- 
tude. The spatial division creates an inner 'core' and outer 
'halo' sample of the prospective candidate targets: the di- 
vision between the core and the halo being arbitrarily de- 
fined as a clustocentric radius of Vp = 30' on the sky, about 
a 4-5 Mpc radius for a typical LARGS cluster. The divi- 
sion in brightness creates a 'bright' and a 'faint' sample. In 
terms of absolute magnitudes, bright galaxies are defined to 
have Rbright < R{M* + 1) and faint galaxies in the inter- 
val R{M* -I- 1) < Rfaint < R{M* + 3). To be consistent 
with P02, M* is chosen to be My = —21.8, this is roughly 
R = 16.5 for a typical LARGS cluster. 

Hence, the final list of spectroscopic targets for this 
cluster consists of four target sets: a bright core of galax- 
ies (rp < 30', R < R{M* + 1)); faint core galaxies (r^ < 30', 
R(M''-|-1) < R< R{M*+3)y, bright halo galaxies (r-p > 30', 
7? < R{AI* + 1)). The fourth possible set of faint halo galax- 
ies was not used because from experience, at these large 
radii the faint sample will be heavily contaminated with field 
galaxies (Figure 1). 

The fiducial stars used by 2dF for guiding are also se- 
lected from the LARGS catalogue. To generate a list of 
suitable fiducial stars the following selection criteria are im- 
posed: CLASS_STAR> 0.9, FWHM < 2.0 and 13 < i? < 15. 
This generates a list of in excess of 30 bright stars. Each 
of these stars are visually inspected to ensure they have no 
diffraction spikes which may cause poor ccntroiding. Their 
astrometric solutions are also compared with the USNO 
Al.O astrometric catalogue (Monet 1996) to ensure that any 
stars with laxge proper motions are removed. 

The allocation of 2dF fibres to our spectroscopic 
targets is done using the CONFIGURE software pack- 
age, supplied with the standard 2dF software (see 
http://www.aao.gov.au/2df/manual/). The spectroscopic 
targets are assigned a priority for the allocation process such 
that n(Bright Gore) > n(Faint Gore) > n(Bright Halo) 
where n(A) is the integer priority given to target set A. Sky 
fibres are allocated and checked to ensure that they are in- 
deed blank sky and not accidentally placed upon an object. 

The fully automated 2dFdr reduction pipeline is used to 
reduce these data. A full description of the package is given 
in Bailey et al. (2001; see also http://www.aao.gov.au/2df/; 
GoUess et al. 2001). 

2.3 Redshift Determination 

Redshift determination is carried out using the XCSAO task 
(Kurtz ct al. 1992) within IRAF and crosscorrelating our 
spectra with numerous template spectra from our private 

§ Galaxies with a flag parameter < 4 are not saturated, merged 
with a secondary source or close to the edge of the CCD. 



collection. The template spectra comprise several stellar 
spectra including synthetic spectra comprising mixes of G- 
and K-type and A- and K-type stellar spectra, a globular 
cluster spectrum and several galactic spectra including spi- 
ral and elliptical galaxies. To determine the reliability of the 
redshift measurements and which of the crosscorrelated red- 
shifts to use we follow the method of Tonry & Davis (1979). 

Tonry & Davis (1979) define a quality value, -Rtdr , as 
the ratio of the height of the fitted peak to the average height 
of peaks in the anti-symmetrical part of the cross correla- 
tion function. The larger the value of Rtdr, the better the 
fit of the template spectrum to the galaxy spectrum. Since 
each 2dF fibre is cross-correlated with many templates, it is 
necessary to determine which of the redshift estimates arc 
correct. The redshift is usually taken to be the value with 
the highest Rtdr- 

To check this, each spectra is de-redshifted to rest 
frame wavelengths and inspected by eye by > 2 authors, 
to check that its emission and absorption features confirm 
the redshift determination (Figure 2). Gonsistent with pre- 
vious investigations (e.g. Golless et al. 2001), we find that 
Rtdr = 3.0 to be a dividing line between good and poor 
redshift determination. We further find that ~ 90 per cent 
of the redshifts found by the above method are verified as 
correct. The remaining fraction have their redshifts manu- 
ally computed. These spectra typically exhibit much noisier 
features, are fainter objects {R ~ M* + 3) and / or are high 
redshift, z > 0.25 and not relevant to this study. Since XC- 
SAO is configured to have an initial guess of a ~ 0.1 for all 
galaxies, these high redshift galaxies often suffer the worst 
determinations as xcSAO attempts to fit the templates to 
incorrect absorption features. Often, a higher initial guess 
{z ~ 0.25) readily solves this problem. 

2.3.1 Internal checks 

We have checked our redshift determinations in two ways: 
firstly by using the emsao task (Mink & Wyatt 1995) to 
generate an independent redshift estimate and secondly by 
using several repeat observations of our clusters. 

EMSAO searches for any emission lines in a given spectra 
by taking the xcSAO determined velocity as an initial guess. 
It then determines a redshift by the fitting of parabolic 
curves to any emission lines found. Since not all of the galax- 
ies have emission lines, this check is only valid for about 20 
per cent of the observations. The median rms of 65kms~^ 
between the independent redshift determination methods is 
very similar to the value of 64kms~^ reported by the 2dF 
Galaxy Redshift Survey (2dFGRS) team (Golless et al. 2001) 
in their repeat observations for their highest quality spectra. 
Further, we find this value docs not depend upon the red- 
shift of the galEixy. Meanwhile, the pair-wise blunder rate, 
defined as those galaxies whose independent redshift deter- 
minations differed by more than 600kms"\ is 2.7 per cent. 
This value is slightly less than the pair-wise blunder rate 
of 3.1 per cent reported by Golless ct al. (2001) for repeat 
observations made by the 2dFGRS team. 

2.4 Success of Star-Galaxy Separation 

For the majority of our clusters, our star-galaxy separation 
techniques presented in POl have been highly successful: the 
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Table 1. Spectroscopic observations undertaken using 2dF. The total number of fibres used per observation is given as N(Fib). This is 
broken down into sky fibres (N(Sky)), spectroscopically confirmed stars (N (Stars)) and galaxies (N(Gal)). 



Cluster Date Grating Seeing Tsxp N(Fib) N(Gal) N(Stars) N(Sky) 

{") (sec) 



Abell 


22 


Jul 08 2000 


600 


V 


1.2-1.6 


5 


X 


1800 


379 


345 


8 


26 


Abell 


1084 


May 07 2002 


600 


V 


1.4-1.5 


5 


X 


1800 


383 


340 


3 


40 






May 08 2002 


600 


V 


1.4-1.6 


5 


X 


1800 


359 


290 


12 


57 


Abell 


1437 


May 04 2002 


600 


V 


1.2-1.4 


5 


X 


1800 


386 


338 


3 


45 






May 05 2002 


600 


V 


1.2-1.3 


5 


X 


1800 


386 


339 


1 


46 


Abell 


1650 


May 07 2002 


600 


V 


1.4-1.5 


5 


X 


1800 


384 


341 


1 


42 






May 08 2002 


600 


V 


1.4-2.0 


5 


X 


1800 


363 


314 


3 


46 


Abell 


1651 


May 04 2002 


600 


V 


1.2-1.4 


5 


X 


1800 


351 


311 


1 


39 






May 05 2002 


600 


V 


1.2-1.3 


5 


X 


1800 


343 


294 


4 


45 






May 08 2002 


600 


V 


1.2-2.0 


5 


X 


1800 


294 


235 


4 


55 






May 08 2002 


600 


V 


1.2-2.0 


5 


X 


1800 


310 


256 


2 


52 


Abell 


1664 


May 06 2002 


600 


V 


1.2-1.5 


5 


X 


1800 


331 


275 


8 


48 






May 06 2002 


600 


V 


1.2-1.5 


5 


X 


1800 


327 


267 


9 


51 


Abell 


2055 


May 25 1998 


600 


V 


1.8-2.2 


4 


X 


1800 


367 


99 


231 


20 






May 28 1998 


600 


V 


1.9-2.5 


5 


X 


1800 


373 


188 


84 


20 


Abell 


2104 


May 26 1998 


600 


V 


2.0-2.4 


5 


X 


1800 


377 


209 


65 


20 






May 27 1998 


600 


V 


1.7-2.2 


4 


X 


1800 


368 


197 


78 


22 


Abell 


2204 


Jul 08 2000 


600 


V 


1.5-1.6 


5 


X 


1800 


379 


280 


65 


34 






Jul 09 2000 


600 


V 


1.7-1.9 


5 


X 


1800 


374 


274 


61 


39 






Jul 09 2000 


600 


V 


1.8-2.0 


5 


X 


1800 


369 


236 


97 


35 


Abell 


3888 


May 26 1998 


600 


V 


1.8-2.3 


5 


X 


1800 


364 


242 


78 


23 






May 27 1998 


600 


V 


2.1-2.4 


5 


X 


1800 


196 


144 


16 


28 






May 06 2002 


600 


V 


1.3-1.5 


5 


X 


1800 


345 


270 


37 


38 


Abell 


3921 


Jul 08 2000 


600 


V 


1.6-1.7 


5 


X 


1800 


366 


331 


3 


32 






Jul 20 2001 


300 


B 


1.7-1.9 


5 


X 


1800 


374 


338 





36 



number of stars amount to only ~few per cent of the total 
number of fibres used (Table 1). There are two exceptions to 
this: the May 1998 observing run and Abell 2204. The May 
1998 observing run used criteria that wore much loss conser- 
vative than those subsequently used (POl; c.f. O'Hely 2000); 
specifically, SExtractor's CLASS_STAR parameter (Berlin & 
Arnouts 1996) was set at a more tolerant level which resulted 
in much higher stellar contamination. Our aim in using such 
tolerant criteria was to try to create a sample with a very 
low bi£is by including more compact objects; akin to the ap- 
proach that Drinkwater and collaborators have successfully 
employed to discover populations of galaxies masquerading 
as stars (e.g. Drinkwater et al. 2003). However, as shown 
in Table 1, this approach proves to be obsorvationally very 
expensive. Subsequent runs therefore adopted the more con- 
servative POl criteria and are accordingly much less contam- 
inated (c.f. Abell 2055 and Abell 1664 for example). From 
Abell 2055, we compute that tiicrc arc 7 compact galaxies 
(of which only 1 is a cluster member; or about 1 per cent of 
the total cluster population) observed with the more toler- 
ant criteria that would not have been observed had we used 
the more conservative criteria. 

The second exception to our success is found in 
Abell 2204 which has a stellar contamination rate of about 
25 per cent (Table 1). Although we use the same criteria as 
POl for selecting the galaxies, we find many more stars in 
it than in other clusters observed during the same run (c.f. 
Abell 22). The reason for this is the line-of-sight proximity of 
Abell 2204 to the Galactic Bulge. Indeed, by applying a cut 
of CLASS-STAR > 0.98 to the photometric parent catalogue 



of Abell 2204, we find that the number of stars grows to over 
eight times the level found in the other parent catalogues at 
faint-end magnitudes {R > 19.0). 

2.5 Line Strength Measurement 

The final stage in the data reduction pipeline is to mea- 
sure various spectral line strengths for each of the galaxies 
observed. The measurement of spectral line strengths is ac- 
complished using an automated program originally written 
by Dr. Lewis Jones. The documentation for the algorithms 
that the program uses are presented in Jones & Worthey 
(1995) and Jones & Couch (1998) and are briefiy summar 
rized here. 

The program calculates the equivalent widths (or 
equally, spectral indices) of up to 30 spectral features by 
finding the ratio of fiux in a 30-40A index passband cen- 
tered on the index to that of a continuum level. To calculate 
the continuum level, the program finds the mean fiux in two 
30-40A sidebands blueward and redward of the index pass- 
band. The local continuum level is then defined by a linear 
interpolation of the mean of these two sidebands. Following 
Jones & Couch (1998), the fiux in each index band is then 
computed as: 

where Fx is the summed flux for a given index, A, Fx,k 
is the fiux at point k for the index and a\^k is the Poisson 
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Figure 2. Examples of combined spectra (left-hand column; created in IRAF and normalized at 4000 Ang) and a small selection of 
individual spectra (right-hand column) that make up the combined ones; split upon spectral classification (see Section 4.3). The spectra 
have been de-redshifted to rest frame wavelengths and several key spectral lines (vertical dotted lines) are indicated. The upper, combined 
a-t-k spectrum is of those galaxies with colours consistent with the CMR. Note that the spectra are not fluxed. 

error bar. The cumulative error bar for a given index, ax, is 
then evaluated as: 

2 _ 1 whore Fc is the flux of the local continuum. Thus the 

"^^ ^1 /ctJ j. resulting values for the measured EW will be negative for 

spectral features that are in emission whilst conversely, pos- 
The equivalent width, EW, for a given index is then: itive values are absorption feature. In this work, all EW 



/ 



(1 - Fx/Fc) 
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values are presented as positive values and are specified as 
being in emission or absorption. 

The detection limit of the measurements of [Oii]A3727A 
is computed by finding the scatter of the EW of the spec- 
tral index about zero and assuming that the median of the 
negative tail of this distribution, 1.5A, is then taken as the 
detection limit estimate. Any modulo value less than this is 
assumed to be noise. 



2. 5. 1 Resolution effect of 300B grating 

Our observations have been made in both main scheduled 
telescope time and service time at the AAT. This has re- 
sulted in two different spectral gratings being used: 600V 
and 300B (see Table 1). The 600V grating has been used 
with the main scheduled telescope time; the 300B grating 
has been used in the service run of Abell 3921 where the 
grating choice was dictated by other service observing pro- 
grammes. 

To assess the effect of the different gratings on subse- 
quent measurements, a selection of bright {R < 19) observa- 
tions made with the 600V gratings are convolved to the reso- 
lution of the 300B observations using the GAUSS task within 
IRAF. The equivalent widths of the H5 and [Oll]A3727A 
spectral line features in the convolved spectra are then re- 
measured. Wc find a median EW offset of 0.46 and -2.33 A 
in and [Oii]A3727 respectively which we herein apply to 
the 300B observations. We emphasize that this offset is only 
relevant for one of our datasots and that the inclusion or 
exclusion of these data from our analysis does not change 
any of our qualitative conclusions. 




16 17 IB 19 30 

R Magnitude 



Figure 3. Tfic completeness for our 2dF programme. The com- 
bined likelifiood (Icft-fiand axis; fiistograms) is the product of a 
given galaxy being botli observed and generating a reliable red- 
shift (-Rtdr > 3.0). The total number of galaxies that this trans- 
lates into is displayed as the curves (right-hand axis). The solid 
curve and histogram refer to those galaxies within 30' (approxi- 
mately 5 Mpc at z ~ 0.1) of the cluster centre whilst the dotted 
ones are for those outside 30'. 



3 COMPLETENESS 

The selection of spectroscopic targets from our photometric 
catalogues is detailed above (Section 2.2) whilst our photo- 
metric completeness limits (e.g. surface brightness consider- 
ations) are dealt with in POl. Briefly, the individual clusters 
do possess different limiting magnitudes in our survey with 
the higher redshift clusters possessing necessarily deeper ob- 
servations. Importantly, the M* + 3 spectroscopic selection 
limit is at least 3(7 above the 5cr photometric completeness 
limit (POl; i.e. all objects are at least 8a detections). 

Moreover, wc do not impose any colour-dependant re- 
quirement upon our selection - the proportions of blue (say 
B — R < 1.6) galaxies to redder ones {B — R > 1.6) in our 
photometric catalogue are statistically the same as the pro- 
portions that are observed with 2dF. Here, therefore, we 
confine ourselves to computing the redshift completeness 
function. 

Since each cluster has had a different number of ob- 
servations made of it with varying degrees of success (i.e. 
successful rcdshifts produced), it is necessary to define a 
completeness function for the spectral sample and take this 
into account when we compare the various cluster popula- 
tions to avoid any possible bias in comparing galaxies at 
small and large radii; comparing our sample with those at 
higher redshifts; or intra-comparing our clusters. 

Two questions have to be answered in order to create 
the completeness function: 



(i) out of the photometric galaxy catalogues, how many 
galaxies were selected for 2dF observation?; 

(ii) out of the galaxies that wore actually observed, how 
many obtained reliable redshift estimates (-Rtdr > 3.0)? 

The completeness function is then simply the product of 
these two probabilities. Figure 3 displays the overall com- 
pleteness for our survey. To generate this plot, we have com- 
bined all of our 2dF observations together and shifted them 
to a common redshift so that their values of M* + 1 is coinci- 
dent. Our coverage of the core, bright samples (radius< 30' 
and magnitudes brighter than M* + 1) is the most complete 
(over 50 per cent of all galaxies brighter than R « 16.5 have 
been observed and produced a reliable redshift; Figure 3), 
followed by the halo bright sample and then the core faint 
sample. Between the individual clusters, the completeness 
levels are comparable - no two clusters differ by more than 
2(7 in their completeness rate at a given magnitude and tar- 
get set. 



4 DYNAMICAL AND SPECTRAL ANALYSIS 

Wo have over 7500 individual spectra, of which about 5100 
are galaxies with good redshift measurements. In this section 
we work out in detail which of these galaxies are members 
of our target clusters. 
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Figure 4. Velocity histograms of the 11 fields examined in this study. In constructing them only galaxies with iixDR > 3.0 have been 
used. The inset panels display an enlargement centred upon the cluster region. 
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4.1 Velocity Structure 

There are a number of ways to qualitatively and quantita- 
tively examine the spectroscopic data in order to determine 
cluster membership. We commence our evaluation of the 
2dF spectroscopy in a qualitative manner. 

We first construct velocity histograms for each of our 
clusters (sec Figure 4). These histograms arc binned into 
1000 kms^^ intervals with an inset panel showing an en- 
largement of a 20000 kms~^ range around the cluster veloc- 
ity, binned into intervals of 400 kms~^. We use these velocity 
histograms to make an initial estimate of the cluster red- 
shifts, indicate possible sub-structure and to identify any 
other obvious structures appearing along the line of sight 
(Appendix B). 

4.2 Cluster Membership via Simple Methods 

We now proceed by defining cluster membership for our 
galaxies. One method to determine cluster membership is 
to use a statistical clipping technique (Colless 1987; C87; 
Zabludoff, Huchra and Geller 1990; ZHG). In their study 
of the uniformity of cluster velocity distributions Yahil and 
Vidal (1977; YV77) defined cluster membership using a 3a" 
clipping technique upon which C87 and ZHG also base their 
methods. YV77 assumes that the velocity distribution of 
the galaxies will follow an underlying Gaussian distribution. 
This is a valid assumption to make if the clusters arc relaxed 
isothermal spheres and the galaxies are massless tracers. The 
relaxed isothermal sphere assumption, however, would be 
invalidated if there is substructure within the clusters (e.g. 
Zabludoff and Zaritsky, 1995). 

Table 2 presents the recession velocities of the clusters 
and their velocity dispersions using ZHG along with other 
global cluster parameters. The values for the errors on the 
velocity dispersion and mean velocity are taken using the 
formula presented in Danese, De Zotti and di Tullio (1980; 
DDD). As with YV77, DDD is vaUd only for Gaussian dis- 
tributions. 

The velocity dispersions for these clusters are in the 
range 700 < az < 1350 km s"^ (Table 2), indicating that 
these clusters are amongst the most massive of their kind at 
these rcdshifts. There is, however, only a very weak (posi- 
tive) correlation between Lx and az- Indeed, a fixed az at 
all Lx is consistent with our data. This is unsurprising given 
the narrow and relatively small range of Lx studied in com- 
parison to other works such as Edge & Stewart (1991) and 
Ortiz-Gil et al. (2004). Objectively, wo can state that within 
the narrow Lx range spanned by the LARCS clusters, the 
scatter in az is ±400 km at a fixed Lx- Ortiz-Gil et al. 
(2004) find that for a much larger sample of 171 REFLEX 
clusters (see Bohringer et al. 2001) an intrinsic scatter of 
az ~ 200 km s~^ and a more significant slope. 

4.3 Cluster membership via a Mass Model 

The statistical clipping techniques used for defining cluster 
membership presented above only makes use of one parame- 
ter: the measured recession velocities. A mass model is capa- 
ble of fully exploiting both spatial and redshift information 
to define the limits of a cluster. The mass model considered 
here is that used by the Canadian Network for Observational 



Cosmology team (CNOC; e.g. Carlberg et al., 1996; Balogh 
ct al., 1999). 

The CNOC mass model derives from a theoretical 
model based upon the assumption that clusters are singular 
isothermal spheres. The technique is described in detail in 
Carlberg, Yee and EUingson (1997; CYE) and briefly sum- 
marized here: Firstly, the difference in velocity. Aw, between 
each galaxy and the mean velocity of the cluster, cz, is com- 
puted. The values of Av are then normalized to the velocity 
dispersion of the cluster, az, and plotted against the pro- 
jected radius away from the centre of the cluster in units of 
^"200^. The mass model of CYE is then used to mark upon 
this plane the 3a and 6a contours (sec especially Table 2 
of CYE and Figure 2 of Carlberg et al. 1997) which are to 
differentiate between cluster galaxies, near-field galaxies and 
field galaxies. 

Using the assumption that a cluster is a singular isother- 
mal sphere, Carlberg et al. (1997) derive the form of r2oo as: 

(3) 5 a, 
= 10^ 

where az is the velocity dispersion of the cluster and 
H{zf = Hi{l + zf{l + noz). 

The results of this analysis are presented in Figure 5. 
The method of ZHG generally agrees well with the CYE 
mass model. Only a small number of galaxies in all of 
the clusters (< 3 per cent) have had their membership re- 
assigned. The majority that are re-assigned now find them- 
selves moved into the category of 'near-field', between the 
3a and 6a contours of the CYE mass model. Furthermore, 
the change in cz and az by using this technique is insignif- 
icant compared to ZHG. Prom herein, therefore, we utilize 
the CYE mass model to define cluster membership. This 
gives us some 1667 cluster galaxies, of which about 1000 arc 
outside r2oo (Table 2). Importantly: we have achieved our 
goal of an average of 150 galaxy members per cluster. 



4.4 Spectroscopic Typing 

Previous work on spectral evolution modeling has shown 
that there are features within a galaxy's spectrum which can 
be readily correlated with their star formation histories (e.g. 
Bruzual & Chariot 1993; Poggianti et al. 1999; Lewis et al. 
2002; Gomez et al. 2003; Balogh et al. 2004). Measurement 
of spectral lines has led some authors (e.g. Couch & Sharpies 
1987; Dressier et al. 1999; Poggianti et al. 1999) to adopt a 
spectroscopic classification scheme based upon the relative 
strengths of certain spectral features. 

Here, the spectral nomenclature of the MORPHs col- 
laboration is adopted (e.g. Dressier et al. 1999; Poggianti et 
al. 1999). This scheme evolved from the work of Couch & 
Sharpies (1987) and Dressier & Gunn (1992). In this scheme, 
galaxies are typed into a particular spectroscopic class on the 
basis of equivalent widths of their H(5A4101 A absorption line 



^ The radius of r200 is defined to be tlic clustoccntric radius at 
which the mean interior density is 200 times tlie critical density 
at the redshift of the cluster, 200pcrit(^)- In terms of Mpc, r200 is 
numerically similar to r^ir which is used by other authors (Girardi 
et al. 1998). 
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Figure 5. The cluster membership technique of CYE appUed to our cluster sample. The solid curve denotes the 3(T contour of the CYE 
mass model; the dashed curve is the dcr contour. An equivalent scale to r/r200 is provided in Mpc along the upper axis for comparison. 
Filled circles are those galaxies defined to be cluster members under ZHG; non-members are open circles. We define ro to be the X-ray 
centre of each cluster. 



© 0000 RAS, MNRAS 000, 000-000 



LARCS III: Spectroscopic Studies 11 



Table 2. Global parameters for the clusters used in this study. We give the right ascension and declination of the X-ray centre of 
each cluster; Lx, the X-ray luminosity in the 0.1-2.4keV passband (Ebeling et al. 1996): the initial clipping used and the results of 
applying the ZHG technique to the LARCS clusters. The errors are calculated from the method of DDD. The final number of galaxies 
with i?TDR > 3.0 present within the cluster is given as N(gal). N(r2oo) is the number of galaxies within r2oo using the CYE cluster 
membership method. Other structures found along the line of sight are noted in Appendix B. 
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against their [Oii]A3727A emission line. In the absence of on- 
going star formation, the strength of the H(5 absorption hne 
is a good indicator of recent star formation within the past 
1-2 Gyr. The presence of [Oll]A3727A in emission is an indi- 
cator of current ongoing star formation in galaxies (Chariot 
& Longhetti 2001). Taken in combination, therefore, these 
spectral line measures can reveal a galaxy's past star forma- 
tion history. The [Oll]A3727A line is chosen as it is readily 
available in the majority of the galaxies. In contrast, Ha (an- 
other potentially useful line to denote current star-forming 
activity; restframe wavelength of 6563A) is unavailable due 
to the spectroscopic range covered. We caution, however, 
that there are drawbacks to using [Oll]A3727A as a star- 
formation indicator; for example, it is sensitive to both dust 
and metallicity of a galaxy. A caveat that also needs mention 
with H(5 is the effect of metallicity (e.g. Thomas, Maraston 
& Korn 2004). H5 is very sensitive to a/ Fe ratio changes at 
super-solar metallicitics and may likely lead to significantly 
younger implied ages for early-type gala^xies. In our analysis, 
we view this effect to be small, as we are only interested in 
comparing similar luminosity systems. 

Our typical lower limit errors on H<5 and [Oii]A3727A 
equivalent widths are 6 and 8 per cent respectively for EW 
> SA, and increasing to 18 and 20 per cent respectively for 
EW ~ 3A. These fi gures are similar to the errors reported 
by Dressier et al. (1999) and hence the application of the 
MORPHs clcissification scheme (Figure 4 of Dressier et al. 
1999) is sensible. There are six different spectral classes in 
the scheme proposed by MORPHs (Poggianti et al. 1999; 
Dressier et al. 1999). This spectroscopic typing scheme is 
described briefly below. 

A graphical outline of the classification scheme is shown 
as Figure 4 in Dressier et al. (1999). The first spectral class 
is simply labelled 'k'. The k class possesses metal absorp- 
tion lines and an absence of - or very weak - emission lines: 
EW([Oii])< 5A and EW(H5)< 3A. They have the same 
properties as the k-type galaxies described by Dressier & 



Gunn (1992): early- type galaxies exhibiting spectra typical 
of K giant stars. The k-|-a class also has metal lines and 
no (or very weak) emission features but possesses moderate 
Balmer absorption: EW([Oll])< 5A and EW(H5) = 3-8 A. 
The a-|-k class possesses stronger Balmer absorption still: 
EW([Oii])< 5 A and EW(H5) > 8A. Together, they replace 
the mixed nomenclature "E-l-A" class (following the sug- 
gestion of Franx, 1993). The letter 'k' here again signifies 
the similarity to K type stars; the 'a' or 'A' correspond- 
ingly signifies similarity to A type stars; i.e. metal lines and 
Balmer absorption present. The k-f-a and a-f-k classes rep- 
resent galaxies that have undergone recent star formation 
that has ceased. Galaxies typed as PSB or HDS by Couch & 
Sharpies (1987) thus now fall into these two classes (Dressier 
et al. 1999). We note that this is only one way of defining 
E-l-A galaxies. Other works (e.g. Zabludoff et al. 1996; Blake 
et al. 2004) use different lines such as H7 to help with clas- 
sification. 

The spectroscopic classes preceded by the letter 'e' 
are galaxies with emission lines present (following Dressier 
& Gunn 1992). Thus the c(c) class has a weak to mod- 
erate H(5 absorption line coupled with a moderate emis- 
sion: EW([Oii]) = 5-40A and EW(H5) < 4A. They are 
well-modeled by an on-going, near constant star forma- 
tion rate (Poggianti et al. 1999). Meanwhile, the e(b) class 
have strong emission indicative of a starbursting galaxy: 
EW[Oii]A3727A > 40A. The e(a) class have a similar 
strength of [Oii]A3727A emission to the e(c) class but ad- 
ditionally possess strong Balmer absorption: EW([Oll]) = 
5-40A and EW(H<5) > 4A. The e(a) and e(b) class can over- 
lap for galaxies possessing EW([Oll]) > 40A and EW(H5) 
> 4A, but this is not an issue for this study as no galgixy 
falls into this classification. See Figure 2 for examples of all 
these spectroscopic types. 

We note that there also exists a small number (< 0.1 
per cent) of further spectroscopic types which cannot be as- 
signed using the two equivalent widths system. One such 
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type are broad- line active galactic nuclei (AGN); these are 
defined by their broad lines and are called e(n) in the 
MORPHs system. Within the 2dF spectroscopy, only one 
broad-line AGN cluster member is found^l although we ac- 
knowledge that there may be many more undetected AGN 
(see Martini et al. 2002). 

Finally, we note that we find a small number of galaxies 
(113) that arc classified as k, k-l-a or a-l-k, yet have significant 
(i.e. EW> 5) H/3 and/or [Oiii]A5007A in emission. 

4.5 Misclassification 

One caveat that should be borne in mind in our spec- 
tral analysis is that LARCS and MORPHs use different 
EW measuring techniques. MORPHs EW measurements 

(Dressier ot al. 1999) arc based upon an interactive Gaus- 
sian line fitting routine which is similar in nature to the 
SPLOT command in IRAF, while we use a passband to de- 
termine the EW of each feature. Using this approach on a 
subset of LARCS data we find that the rms difference in 
measured EW from LARCS to MORPHs is « 0.6A for both 
[On]A3727A and H5. 

The existence of such a difference will have a blurring 
effect on the boundaries between the various spectroscopic 
types. We have investigated how much of an effect this is 
by perturbing the H5 and [Oll]A3727A EW measurements 
100 times in a Monte-Caxlo fashion. We then recompute the 
spectral cleissification using these new values. The majority 
of spectra that get re-typed arc k galaxies scattering into 
k-|-a types, with a smaller number scattering from e(c) to 
e(a) types and k-|-a to k types. We note that at all classifi- 
cation boundaries, the re-typing rate is not greater than 5 
per cent. 

Whilst misclassification could be an issue, we find that 
there are real differences in the typical spectral energy dis- 
tributions between the average k-|-a and k types. In Figure 2, 
we display combined spectra for each class. These average 
spectra show that the EW of the H5 line increases from k 
(1.3 A) through k+a (3.9 A) to a+k (8.9 A). The Balmer 
H7 absorption line also strengthens in our combined spectra 
(Figure 2) from k (0.9 A) through k-|-a (2.0 A) to a-|-k (4.1 
A) types, as expected if these galaxies contain progressively 
larger populations of A- and F-type stars. 



5 ANALYSIS, RESULTS AND DISCUSSION 
5.1 Colour-Magnitude Relation 

The more luminous early-type galaxies (ellipticals and SOs) 
within local clusters arc known to exhibit systematically red- 
der integrated colours than the late-types (Visvanathan & 

II To find one cluster AGN Is not unexpected. Dressier, Tliomp- 

son, & Shcctman (DTS; 1985) report that for a sample of 1268 
galaxies over 14 clusters at z ~ 0.04, 1 per cent (i.e. 12) of all 
cluster galaxies possess active nuclei. Within the LARCS obser- 
vations, therefore, one expects to find ~ 4 AGN. The absolute 
magnitude limit probed by DTS is brighter than LARCS, but 
the number of expected AGN should not increase by extension 
to the fainter LARCS absolute magnitude limit (Huchra & Burg 
1992). 



Sandage 1977). These galaxies also demonstrate a tight cor- 
relation between their colours and magnitudes: the colour- 
magnitude relationship (CMR; e.g. Bower, Lucey and El- 
lis 1992). For sub-L* galaxies, Kodama & Arimoto (1997) 
demonstrate that the slope of the colour- magnitude relation 
(assuming that the stellar populations are formed in a sin- 
gle event) could be due to the mean stellar metallicity (as 
suggested by Dressier 1984), with the scatter around the re- 
lation being due to age effects (Kodama et al. 1999). Fainter 
galaxies on the CMR may have been constructed from the 
fading or disk stripping of later-type galax;ies, thus leaving 
a bulge with a small amount of younger stars present. 

Following earlier studies by Abraham et al. (1996) and 
Terlevich et al. (2001), P02 and Wake et al. (2005) sug- 
gest that the CMR of clusters evolves bluewards with in- 
creasing clustocentric radius and decreasing local galaxy 
density. Their investigations, however, are made in the ab- 
sence of spectroscopy by making a statistical correction 
to remove the contaminating field population. By employ- 
ing spectroscopic membership, colour-magnitude diagrams 
are now constructed for spectroscopically-confirmed cluster 
members. 

We use a biweight method to fit the CMR in an identical 
manner to that used by P02 (see also Beers, Flynn & Geb- 
hardt 1990; Press et al. 1992). The spectroscopically-coded 
colour-magnitude diagrams for our eleven clusters whose 
galaxies are within a projected radius of < 2 Mpc of the 
cluster core are presented in Figure 6 with the biweighted 
fits to their OMRs. 

To proceed, we now exploit the homogeneity of our clus- 
ters and combine them to make a composite. This is per- 
formed in the same manner as P02. Briefly, we have to trans- 
form the colours, magnitudes and positions of the galaxies 
onto a common scale. The apparent magnitudes are evolved 
to a median redshift of z = 0.12 using the relative magni- 
tudes of L* galaxies. The colours arc dealt with by reducing 
the observed CMRs to provide colour information relative 
to the fitted relation and then transformed to z = 0.12. 
Positions are moved onto a common scale by scaling to a 
fixed metric size. The relative merits of these transforma- 
tions are discussed in detail by P02 and Wake (2003). Here 
it is sufficient to note that our clusters are a homogeneously 
selected sample with similar X-ray luminosities (i.e. clus- 
ter masses), therefore normalizing to a fixed metric size or 
to r2oo (Table 2) makes no significant differences to the re- 
sults presented (cf. Wake et al. 2005 whose sample spans a 
factor of over 100 in X-ray luminosity). Finally, we create 
colour histograms of the composite cluster down to a lim- 
iting magnitude of AIv = —20 to search for any variation. 
Since our composite spectroscopic sample is much smaller 
than the full photometric one used in P02, we restrict our 
radial bins to be 2 Mpc annuli. To better identify these his- 
togram's peak values and quantify the radial blueing, we fit 
them with Gaussians using a minimization method. More 
extensive details of our fitting routine can be found in P02. 
These Gaussians are tabulated in Table 3. The peak colour 
of the CMR in the composite cluster evolves blueward at a 
rate of d{B — R) /drp = —0.011 ± 0.003 which is a rate some 
2.2cr smaller than found from the photometric method of 
P02: d{B-R)/drp = -0.022 ±0.004. Therefore those galax- 
ies in the cluster outskirts may have luminosity-weighted 
mean stellar populations up to 3 Gyr younger than those 
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Figure 6. Colour— magnitude diagrams for the 11 clusters used in this work. Cluster members within 2 Mpc of the cluster centre are 
plotted. The biweight fit to these distributions is shown as the solid line; the parallel flanking dashed lines show the la uncertainty in 
the colour. The vertical dotted lines denote My = —20 which is used to define the sample for the biweight fit. The lower right corner 
displays our field sample for comparison. This is constructed from non-cluster members and evolved to z ~ 0.1; the CMR lines are to 
illustrate where the CMR would be (we note that others do report the existence of a CMR for field ellipticals; e.g. see Bell et al. 2004) 
at this redshift; see Figure 2 of P02. The galaxies are coded upon the spectroscopic types: solid (red) circles are k types; open (magenta) 
triangles are k+a and a+k types; starred (blue) points are e types. 



in the cluster centre assuming d{B — R)/dt = 0.03 mag per 
Gyr (Kodama & Arimoto 1997; see also Moran et al. 2005). 
One can also interpret Table 3 as showing that most of the 
radial blueing occurs within the central ~ 3-4 Mpc of the 
cluster centre. This would make the rate much steeper and 
inline with the radial colour changes noted by Gerken et al. 
(2005). It is also possible that at these large radii from the 
cluster centre, we may have a mixed galaxy population of 
those that have already been through the cluster centre and 
those that are infalling for the first time (e.g. Rines et al. 



2005); we will explore this scenario in Section 5 in greater 
detail. 

To examine this trend in a somewhat more general way 
we use local galaxy density, E, instead of clustocentric ra- 
dius. The local galaxy density is estimated from the photo- 
metric catalogue by finding the surface area on the sky that 
is occupied by a given galaxy and its 10 nearest neighbours 
dowrn to Mv = —20. As with the photometric method, this 
value will be overestimated due to background galaxy con- 
tamination, therefore we correct these values by subtracting 



© 0000 RAS, MNRAS 000, 000-000 



14 K.A. Pimbblet et al. 



1.5 



o 







1 - 



D 0.5 



0.8 - 



^0.6 
2 



X 

(0 

■S 0.2 



« 



' ' ' I ' ' ' ' I 

log(S) > 2.0 



T — n — I — I — I — I — I — I I I — I — n 
1.5 < log(£) < 2.0 



T — I — I — I — I — I — n — I — I — TT 
1.0 < log(S) < 1.5 



0.5 




1.5 

(B-R) 



2.5 



Figure 7. Colour histograms of the composite cluster split on local galaxy density, S. The best fitting Gaussians for each bin are 
displayed as the dotted curves. The peak of the richest bin is noted on each panel by the vertical solid line whilst the peaJc of each 
individual bin is shown as the dashed vertical line. A blueward shift in the peak of the OMR is seen with decreasing local galaxy density 
at a rate of d{B - R)/dlog{ll) = -0.062 ± 0.009. 



off a constant density computed from the median local den- 
sity of a field sample (P02). Since the range in S is large, 
we work in logarithmic values and generate six bins covering 
three orders of magnitude in /op(E). 

For each of these bins, we create a colour histogram from 
the composite spectroscopic cluster sample and fit Gaussians 
to them using our minimization method. The result of 
this analysis is presented in Figure 7 and tabulated in Ta- 
ble 3. Here, the peak colour of the CMR evolves blucwards 
with local galaxy density at a rate of d{B — R) / dlogiT,) = 
—0.062 ± 0.009 which is, again, accompanied by a broad- 
ening of the CMR peak. Like the radial blueward trend, 
the trend with local galaxy density is smaller than the 
one found by the photometric method in P02 by ~ Icr: 
d{B - R)/dlog{T) = -0.076 ± 0.009. 

From Figure 7 and Table 3, we are also able to examine 
how the width of the CMR changes with environment. The 
Gaussians fitted to these data are used to compute upeak, the 
CMR's width. The width increases by some '-^0.1 across the 
range studied, consistent with our previous findings (P02). 
We note, however, that even though we are using Gaussians 
to parametrize the CMR's, there is no physical reason why 
the CMR should be well fit by a such a curve. 

To better parametrize the changes in the shape of the 
CMR distribution we look at the trends in the quartile 
points with radius. The reddest galaxies on the CMR, evalu- 
ated by computing the SO**" percentile of the colour distribu- 
tion, possess a near constant {B — R) colour over the range 



Table 3. Peak colour and full width of the variation of the com- 
posite CMR with both clustocentric radius and local galaxy den- 
sity as derived from fitting Gaussians to the colour distributions 
(see also Figure 7). 



Sample Peak (B - R)my=-21.S crpeak 



Radius (Mpc) 

0-2 1.86 ± 0.01 0.12 ± 0.02 

2-4 1.78 ± 0.03 0.14 ± 0.03 

4-6 1.77 ± 0.04 0.16 ± 0.05 

6-8 1.77 ± 0.07 0.21 ± 0.06 

logio (Local Galaxy Density) 

> 2.0 1.87 ± 0.01 0.12 ± 0.02 

1.5-2.0 1.83 ± 0.02 0.14 ± 0.03 

1.0-1.5 1.84 ± 0.04 0.15 ± 0.05 

0.5-1.0 1.80 ± 0.05 0.15 ± 0.06 

-0.5-0.5 1.78 it 0.07 0.16 ± 0.06 

< -0.5 1.77 ± 0.09 0.18 ± 0.08 



of projected radius and local galaxy density probed in this 
study. Conversely, the bluer members, examined using the 
70*'' percentile, exhibit a strong blueward shift. This result 
is very similar to the P02 result (see Figure 7 of P02; also 
see Figure 13 of Wake et al. 2005). 
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Figure 8. Colour-magnitude diagram for the composite cluster coded upon spectral type. The vertical dashed line denotes the fiducial 
magnitude cut My = —20 that we use. The horizontal solid line is the position of the composite cluster's best fit CMR (note that the 
slope of the CMR has been removed in producing the composite cluster). 



One caveat in this analysis is that the best sampled sys- 
tems such as Abell 3888 possess a greater weight in the com- 
posite cluster than the poorest sampled (e.g. Abell 1664). To 
explore this issue in more depth, we normalize all the clus- 
ters to have the same sampling and re-construct our com- 
posite cluster. Repeating the above experiments results in a 
change of rate in d{B — R) /drp and d{B — R)/dlog{T,) by no 
more than 1.5a - i.e. no significant change. Therefore, whilst 
there are cluster-to-cluster variations in this sample, we em- 
phasize that the LARCS sample is an homogeneously X-ray 
selected one, and such variations are therefore minimized. 



5.2 Populations on the CM plane 

With the peak colour of the CMR progressing blueward with 
radius and local galaxy density, we may expect that the 
mean colour of the passive k-type galaxies on the CMR will 
also reflect these changes. 

Figure 8 shows the colour magnitude diagrams of the 
different spectroscopic types in the composite cluster. His- 
tograms of the colour distribution are shown in Figure 9 
down to Mv ~ —20. These colour distributions show an 
number of interesting features. In the k-type class, there 



are a number of apparently non-starforming galaxies with 
colours bluer than the CMR and a small number of redder 
ones. The bulk of the k-types, however, do possess colours 
consistent with the CMR; ~ 85 per cent of the population 
lies in the range 1.6 < (B—R) < 2.0. Although there appears 
a large scatter in the k-type galaxies, it is no more than the 
scatter from the individual clusters (Figure 6). The post- 
starburst types (k-l-a and a-l-k) have a much flatter {B — R) 
distribution. Whilst most lie blueward of the CMR, approx- 
imately 40 per cent of them have colours consistent with 
the CMR. We note, however, that some of the poststarburst 
types with colours consistent with the CMR are close to 
the boundary with k types (i.e. EW(H5)« SA) and may 
therefore have been misclassified (see above). Also, the k-l-a 
colour distribution appears to have a much larger scatter 
than the a-|-k one. However, with the small number (20) of 
a-|-k galaxies available to us (Figure 6), any firm conclusion 
about colour differences between the a-|-k and k-l-a popula- 
tions is likely premature. The emission line types (e(a) and 
e(c)) are qualitatively different distributions, exhibiting sig- 
nificant fractions of blue {{B — R) < 1.6) galaxies (55 and 
90 per cent respectively). Indeed, the e(a) class is by far the 
bluest and faintest class of galaxy reflecting its characteris- 



© 0000 RAS, MNRAS 000, 000-000 



16 K. A. Pimbblet et al. 



tic vigorous on-going star-formation, whilst the e(c) types 
have a colour distribution similar to the k-f a types. 

On the issue of galaxy colours, we noted that Figure 8 
shows a small yet significant population of blue passive (k 
type) galaxies {B — R < 1.6) and very red galaxies {B — R> 
2.0; i.e. redder than the CMR ridge). To investigate the stel- 
lar populations of these galaxies, we combine together mul- 
tiple spectra using the IRAF task scombine to improve the 
signal to noise ratio instead of using individual EW mea- 
surements from each galaxy (Figure 10). In the combined 
spectra, each individual spectra is weighted according to our 
completeness function in order to remove luminosity differ- 
ences. 

We find that the shape of the combined blue k-type 
spectra (bottom of Figure 10) is much more like a typical 
spiral galaxy than for the CMR sample. But compared to 
those galaxies on the CMR, the blue k-type galaxies have 
a small amount of [Oll]A3727A present coupled with obvi- 
ous H/3 and [Oiii]A5007A crriissiori lines. The EW of this 
[Oll]A3727A is 2.6A and hence on average these galaxies 
have 1.5 to 2 a [Oii]A3727A emission individually which 
in the MORPHs spectroscopic classification system ensures 
that they are typed as k's. Meanwhile, the EW of the H/3 
and [Oiii]A5007A Unes are 2.6 and 2.8 A respectively. We 
suggest that these passive blue galaxies have suffered some 
recent star formation which has perturbed their luminosity- 
weighted broadband colours. If their modest star formation 
rate ceases, then they arc likely to quickly evolve on to the 
CMR (see Figures 6 and 7 of Kodama & Bower 2001). 

The very red galaxies (top spectra of Figure 10) are 
similar to those on the CMR,; lacking obvious [Oll]A3727A 
emission combined with small amounts of US. Therefore, we 
tentatively suggest that these very red cluster galaxies have 
colours that are rather dust affected or are super-metal-rich. 
If they are dust affected, then it would not only affect their 
colours, but also produce an absence of UV emission lines. 

Although we have no detailed morphological classifica- 
tions for our sample (Pimbblet 2001), we can use the dis- 
tribution of our galaxies on a plane of concentration in- 
dex (CI) and maximum surface brightness (^imax) to look 
for changes in the morphological mix within the compos- 
ite cluster (POl; Abraham et al. 1994). On the CI-yUMAx 
plane, early-type galaxies populate the high-concentration 
and high surface brightness region, with later-types typi- 
cally having lower concentrations and fainter maximum sur- 
face brightnesses (POl). Comparing the distribution in CI- 
A*MAx of galaxies, we find that the very red sample is sta- 
tistically the same as those galaxies drawn from the CMR. 
However, the blue k-type sample is composed of many more 
low /Umax plus low CI galaxies than the CMR sample is. Ap- 
plying a two dimensional Kolmogorov-Smirnov (K-S) test 
(Fasano & Franceschini 1987) to the CI-/Umax distribution 
we find that the blue k-types are unlikely to have been drawn 
from the same parent population as the CMR sample at a 
98 per cent confidence level. Since late-type galaxies in the 
field typically have much stronger [Oll]A3727A (Jansen et al. 
2000) for their observed continuum colours, the blue k types 
are not normal disc systems. We suggest that these galaxies 
may therefore be anaemic late- types (i.e. similar to late-type 
galaxies except they possess low or insignificant amounts of 
emission) analogous to those observed at higher redshifts by 
Poggianti et al. (1999) and Balogh et al. (2002). At most. 
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Figure 9. Colour distribution of the different spectroscopic types 
contained within the composite cluster from Figure 8 with a mag- 
nitude cut-ofT of My = —20. 

they comprise 2 per cent of the whole cluster population 
at ~ 0.1 which is somewhat lower {2a significance) than 

the 6.5 per cent reported by Balogh et al. (2002) across 10 
z 0.25 clusters for their anaemic late-type population. 

5.3 Spectral populations in individual clusters 

We present the breakdown by spectral class of the clus- 
ter members and a field sample** in Table 4. The major- 
ity of the clusters have a dominant k component, typically 

accounting for half to three-quarters of the entire cluster 
population to L* + 2. There are a significant number of 
post-starburst galaxies (k-|-a and a-|-k), accounting for a fur- 
ther ~10 to 20 per cent. This is in contrast with the recent 
study of the similarly X-ray luminous local Coma cluster by 
Poggianti et al. (2004) who find no post-starburst galaxies 
at such bright magnitudes. Moreover, unlike Poggianti et 
al. (2004), our post-starburst galaxies are not preferentially 
situated anywhere in particular in the cluster, regardless of 
their colour. The actively starforming^t fraction tends to be 

** Our field sample is generated by combining those galaxies that 
are non-cluster members and not contained in secondary struc- 
tures (Appendix B) in the redshift range 0.07 < z < 0.15. 
tt Here, we use the term actively starforming to refer to galaxies 
that possess certain emission line properties. We intentionally do 
not state whether these are star forming or AGN or both as we 
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Table 4. Spectroscopic classifications of cluster galaxies using the MORPHs system. Note that we find no e(a,b) types in our clusters 
while only A1664 and A2055 have any e(b) type galaxies. The numbers are all percentages and completeness function weighted (Figure 3). 
The final rows are for our composite cluster and field sample: constructed from the individual cluster observations and weighted by the 
completeness function. The column headed 'Cluster Morphology' gives a coarse indication of the general state of the cluster from work 
presented in P02 and from visually inspecting archival ROSAT X-ray imaging of our clusters. 



Cluster 


n(gal) 






Percentage 






Cluster 






k 


k+a 


a+k 


e(c) 


e(a) 


e(b) 


Morphology 


Abell 22 


67 


54 ± 5 


15 ± 3 


± 1 


18 ± 3 


12 ± 2 


nil 


Regular to intermediate cluster 


Abell 1084 


122 


59 ± 4 


19 ± 2 


1 ± 1 


12 ± 2 


9 ± 2 


nil 


Irregular cluster 


Abell 1437 


224 


67 ± 3 


8 ± 1 


4 ± 1 


11 ± 1 


10 ± 1 


nil 


Regular cluster 


Abell 1650 


208 


67 ± 5 


9 ± 3 


± 1 


14 ± 3 


10 ± 2 


nil 


Regular cluster 


Abell 1651 


208 


75 ± 7 


14 ± 2 


3 ± 1 


6 ± 1 


2 ± 1 


nil 


Regular cluster 


Abell 1664 


123 


51 ± 4 


10 ± 2 


2 ± 1 


20 ± 2 


15 ± 2 


1 ± 1 


Highly irregular cluster 


Abell 2055 


107 


77 ± 7 


15 ± 3 


± 1 


7 ± 2 


1 ± 1 


± 


Regular cluster 


Abell 2104 


156 


47 ± 4 


25 ± 3 


7 ± 2 


15 ± 2 


4 ± 2 


nil 


Intermediate cluster 


Abell 2204 


125 


66 ± 7 


11 ± 3 


5 ± 3 


16 ± 4 


1 ± 1 


nil 


Intermediate cluster 


Abell 3888 


201 


70 ± 5 


12 ± 2 


3 ± 2 


10 ± 2 


7 ± 2 


nil 


Intermediate cluster (see Girardi et al. 1997) 


Abell 3921 


126 


56 ± 5 


20 ± 4 


1 ± 1 


14 ± 3 


6 ± 2 


nil 


Undergoing merger event (Ferrari et al. 2005) 


Composite 


1667 


63 ± 2 


14 ± 1 


2 ± 


13 ± 1 


7 ± 1 


± 


N/A 


Cluster < r200 


692 


71 ± 3 


10 ± 1 


1 ± 


10 ± 1 


3 ± 1 


nil 


N/A 


Cluster > r2oo 


975 


59 ± 3 


16 ± 1 


1 ± 


16 ± 1 


9 ± 1 


± 


N/A 


Field 


853 


48 ± 3 


20 ± 1 


3 ± 1 


18 ± 2 


6 ± 1 


± 


N/A 
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Figure 10. Restframe combined k type spectra (non-fluxed) for 
galaxies in the composite cluster with {B — R) < 1.6 (bottom), 
1.6 < [B ~ R) < 2.0 (centre) and {B - R) > 2.0 (top). Key 
spectral features are marked (top and dotted lines). 

marginally greater than the poststarburst fraction, amount- 
ing to ~15 to 25 per cent of most of our clusters. 

In contrast to the MORPHs clusters (Table 4 of Pog- 



cannot distinguish between these possibilities (see Martini et al. 
2002). 



gianti et al. 1999), our cluster core regions are much less ac- 
tive than their z ~ 0.4 cousins; i.e. higher passsivc fractions 
and lower emission fractions in general. To better compare 
directly to the MORPHs clusters, we note that the Poggianti 
ct al. (1999) sample is limited to My = —19; roughly equal 
to M* + 1.5 (Dressier et al. 1997) and that their spatial ex- 
tent is roughly 2 Mpc at their median redshift - about or 
just below 7-200 for our clusters. By imposing the same limits 
on our composite cluster sample, the fractions presented for 
the rp < r20() sample in Table 4 do not change by more than 
2a in a given spectral type. We tentatively suggest therefore, 
that the e(b) class has vanished to all intents and purposes at 
z ^ 0.1 whereas Poggianti et al. (1999) find about 5 per cent 
of cluster members are e(b) types at « ~ 0.4. This is only 
tentative as 5 per cent differences between our sample and 
Poggianti et al. (1999) are of marginal significance - about 
2.5a. Further, from z ~ 0.4 to z ~ 0.1, the e(c) fraction 
remains approximately the same; the e(a) fraction falls by 
about 4 per cent (not significant); the combined k-|-a Sz a.+k 
faction falls by 5 per cent (not significant); while the k frac- 
tion increases by some 15 per cent (> 3cr significance). The 
k types may be the endpoints of active galaxies once their 
activity fades and hence it is not surprising that evolution 
is most marked in this subsample compared to the different 
flavours of active galaxies. We also note that MORPHs is a 
mixed cluster sample spanning a range in richness and are 
not all X-ray luminous. 

Although our sample has uniform X-ray selection, there 
also appears to be real variation between our clusters. For 
example, our most irregular cluster, A1664, is marked out 
by a very small ratio of active to passive types (e(a)-l-e(c) 
: k = 1 : 1.5) whereas A2055 and A1651 both have a much 
larger ratio (~ 1 : 9.5). Such diff'crcnces in passive : ac- 
tive ratios can not be wholly accounted for by differences in 
global cluster parameters such as X-ray luminosity or veloc- 
ity dispersion (Table 2). However, the coarse morphologies 
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Figure 11. Fraction of active galaxies (i.e. 1 - fraction of k types) 
in the clusters (Table 4) versus the of a Gaussian fit to their 
velocity distributions Figure 4. The solid line is a fit to the points 
and is 3o" away from a slope of zero. Due to this mild correlation 
between the two variables, we suggest that the clusters with the 
least Gaussian velocity distributions also have the lowest fraction 
of k type galaxies, or equally, are the most active clusters. 



of the clusters (see Table 4; our coarse morphologies are de- 
rived from visually inspecting ROSAT X-ray maps and two 
dimensional galaxy distribution from P02) show that the 
more irregular systems (e.g. A1664 and A3921) tend to dis- 
play a higher incidence of active galaxies than their more 
regular counterparts (e.g. A1437). We note that A1664 is 
likely in the early phases of merging and the enhanced ac- 
tive fraction could be a result of this. While, A3921 is in 
the middle phase of merging (Ferrari et al. 2005) and is less 
active, but displays many more post-staxburst types. More- 
over, the clusters with the smallest fractions of k type galax- 
ies also appear to be the ones with the most non-Gaussian 
appearing velocity distributions (Figure 4). To test this, we 
fit Gaussians to the N{z) distributions shown in the inset 
panels of Figure 4, allowing the normalization, mean veloc- 
ity and a of the fits to all be free parameters, and find the 
of each fit. In Figure 11 we show the relationship between 
the and the fraction of active galaxies in each cluster. 
We find a correlation between these variables, at a 3cr level. 
Therefore, the clusters that possess the least Gaussian ve- 
locity dispersions tend to be the ones that either are actively 
starforming or have recently formed stars. 

We conclude that even in an homogeneously X-ray se- 
lected sample of galaxy clusters, individual clusters can show 
real and significant variation from cluster to cluster, prob- 
ably reflecting their recent infall histories. The only corre- 
lation with fractions of different spectral types is the mild 
anti-correlation between Gaussian velocity distributions and 
k type fraction. We suggest that these variations between 
the clusters are likely driven by local processes as there is 



no clear correlation with other global parameters such as Lx 
or az with the highest / lowest fractions of post-starburst 
or active galaxies. 

5.4 Environmental Dependence of Line Strengths 

The trends in the characteristic colour of the CMR we see 
as a function of radius or local galaxy density should be 
mirrored by changes in the spectroscopic properties of these 
galaxies. In particular, if the variation in the {B — R) colour 
comes about from its sensitivity to the age of the stellar 
population (e.g. Searle et al. 1973), then we would expect 
to see increasing signs of youth in the spectra of galaxies on 
the CMR in the lower density outskirts of the clusters. The 
Balmer absorption line, HJ, provides a useful indicator of 
mean stellar age in non-stax forming systems (although we 
note that does not trace only age; Thomas et al. 2004). 
Indeed, Terlevich et al. (1999) demonstrate the relationship 
between deviations in the {U — V) colour of galaxies away 
from the CMR and the equivalent deviations in the strength 
of Balmer absorption lines from observations of early-type 
galaxies in the Coma cluster. Based on a CMR zeropoint 
shift of A{B — R) ^ 0.07 for typical galaxies out to 4Mpc 
(Table 3), we estimate (using data from Terlevich et al. 1999; 
see also P02) that this population should show an enhanced 
Balmer absorption of ^ 2.8A equivalent width. 

Using only the non-starforming population of galax- 
ies (i.e. the emission line types of galaxies with 
EW[Oii]A3727A > 5.0 are excluded) that he on the CMR, 
a detection of an increase in the mean strength of the US 
absorption line by ~ 2.2A from the cluster core region to 
the outskirts of the cluster is confirmed at > 3.5a (Table 5). 
However, this level of change in H5 could have been expected 
by 4 Mpc rather than in the outskirts (i.e. by 8 Mpc). We 
note that most of the colour change occurs within the cen- 
tral ~ 4 Mpc and further colour changes with radius are 
minimal. Therefore whilst the prediction remains valid at 8 
Mpc; at 4 Mpc, we are at least 2.5a away from the prediction 
of Terlevich et al. (1999). This may be related to the fact 
that Terlevich et al. (1999) study the Coma cluster which 
is at irmch lower redshift than our sample. Fitting a line 
to the points presented in Table 5 weighted by their errors 
yields a rate of dHS/drp = 0.35 ± 0.06. The trend with local 
galaxy density, however, is more noisy. The highest density 
bin (log(S) > 2.0) has a large error due to the comparatively 
small number of galaxies present. Weighting by the errors 
in Table 5 also gives a rate of dR5 / dlogiji) = 0.66 ± 0.15 
with decreasing density. To ascertain whether radius or local 
galaxy density is more important, we compute the relation- 
ship between the two variables. Then, for given radii (or 
local galaxy densities), we find the corresponding density 
(or radius) and compare which of the two rates shows the 
stronger trend. In most cases, local galaxy density appears 
to be more important (stronger) than radius from the cluster 
centre, although the differences at large radii (> 3 Mpc) be- 
come increasingly less significant. Therefore, we tentatively 
suggest that local galaxy density is a more fundamental pa- 
rameter than radius (in contrast, see Whitmore et al. 1993). 
We also note that we have included the more irregular clus- 
ters in this analysis (Table 4); our result does not change 
(i.e. local galaxy density being more important) even if we 
exclude these clusters from the composite. 
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Table 5. Increase of the mean strength of the H5 absorption line 
in passive galaxies with colours consistent with the CMR (within 
±1(t) as a function of both cluster radius and local galaxy density 
in the composite cluster. 



Sample Mean EW(H5) (A) 



Radius (Mpc) 

0-2 0.39 ± 0.16 

2-4 1.06 ± 0.28 

4-6 1.24 ± 0.36 

6-8 2.56 ± 0.69 

log (Local Galaxy Density) 

> 2.0 0.83 ± 0.40 

1.5-2.0 0.36 ± 0.25 

1.0-1.5 1.05 ± 0.32 

0.5-1.0 0.51 ± 0.41 

-0.5-0.5 1.07 ± 0.40 

< -0.5 1.95 ± 0.46 



Gomez et al. (2003) and Lewis et al. (2002; see also 

Balogh ot al. 2004) report that the starformation rate in- 
creases significantly at at large distances from the cluster 
core (and similarly, for small values of local galaxy density) . 
For EW[Oii]A3727A, Gomez et al. (2003) demonstrate that 
the median value starts to differ significantly from the field 
value at about 3-4 virial radii irv) from the cluster core. 
This critical radius corresponds to roughly « 4 Mpc (or 
~ few galaxies per Mpc^) for our survey and assumed cos- 
mology. 

In Figure 12, we plot the median EW [Oii]A3727A for all 
our cluster galaxies brighter than M* + 1 (to match Gomez 
et al. 2003) as a function of rp and E. For comparison, we 
have also constructed a, z Q.l 'field' sample. This is done 
by taking all galaxies in the range 0.07 < z < 0.16 (to 
match the redshift of the LARCS clusters) and then ex- 
cluding all cluster members and other significant structure 
(e.g. the wall in Abell 22; Pimbblet, Edge & Couch 2005; 
Abell 1079 which is in the periphery of Abell 1084; Pimbblet 
& Drinkwater 2004; P02). From this sample, we compute the 
75th and 25th quartiles of the EW [Oii]A3727A distribution 
to compare directly with the cluster sample; see Figure 12. 
The median EW [Oii]A3727A increases from the cluster core 
to the outskirts of the composite cluster, with a 'break' at 
around 3-4 Mpc. The changes in the median shown in Fig- 
ure 12 are due to a decline of the upper quartile. The break 
at 3-4 Mpc corresponds well with the characteristic radius 
found by Gomez et al. (2003). In order to test where the 
cluster population becomes significantly different from the 
field sample, we employ a KS test on the two distributions. 
We find that the cluster EW [On]A3727A distribution starts 
to differ significantly (> Ict; 68 per cent level) from the field 
distribution at below Vp < 3.8 Mpc (Figure 12). Such a ra- 
dius is very comparable to the Gomez et al. (2003) result 
(rp « 4 Mpc; vertical dotted line in Figure 12). In terms of 
local galaxy density, this radius corresponds to « 1.3 ± 0.2 
galaxies per square Mpc (c.f ~ 1 galaxy per square Mpc 
found by Gomez et al. 2003). 



Table 6. Dispersion in Av / cr^, relative to the k types, and 
average distance from the cluster centre per spectral type from 
Figures 14 and 15. The star-forming types tend to reside at the 
outskirts of the cluster and possess a much larger velocity disper- 
sion. 



Type 


a (Av / (Tz) 


radius / r200 


k 


1.00 


1.01 ±0.06 


k-l-a 


1.61 


1.49 ±0.17 


a-l-k 


3.41 


1.73 ±0.23 


e(c) 


3.20 


1.76 ±0.21 


e(a) 


5.85 


2.09 ±0.16 



5.5 Luminosity functions and phase-space 

diagrams 

To better understand the relationship between the different 
spectral classes in clusters and their evolution from higher 
redshift counterparts (i.e. Dressier et al. 1999; Poggianti et 
al. 1999), we now examine their absolute luminosity distri- 
butions (Figure 13) and orbital characteristics (Figure 14). 
The absolute magnitude distributions for k, k+a and a+k 
types are brighter than for the e(c) and particularly for the 
e(a) types (similar to the z ~ 0.4 sample from Dressier et al. 
1999). However, Dressier et al. (1999) report that at 2: ~ 0.4, 
the distributions of k, k+a and a+k are indistinguishable, 
but in our sample at « ~ 0.1, the k+a and a+k luminosity 
function is much flatter than for the k types (Figure 13). 
This may be due to mixing bright k types with fainter, real 
k+a and a+k types as a result of scattering across the spec- 
tral type boundaries (Section 4.5). 

To examine how the different spectroscopic types are 
distributed within the composite cluster, we construct both 
velocity-radius plots of the different spectral types (Fig- 
ure 14) and a cumulative radial distribution, derived for all 
galaxies brighter than My = —20 (Figure 15). The pas- 
sive k-type galaxies dominate the core regions of the cluster 
(r-p < r200) whilst the other types occur more frequently 
beyond r2oo and are generally of higher velocity dispersion. 
To quantify this, we compute a normalized velocity disper- 
sion, cr(Av / ffz), and mean radius of each of the different 
spectral types (Table 6). This demonstrates that the active 
types are much more likely to reside in the outskirts of the 
cluster (similar to the finding of Gerken et al. 2005 that 
the fraction of emission line galaxies increases strongly be- 
yond 2 Virial radii) and possess a much higher velocity dis- 
persion range than the poststarburst or the passive types. 
This is confirmed in Figure 15, where the passive k-type 
galaxies are distinctly seen to be the most centrally concen- 
trated whilst the emission line galaxies preferentially avoid 
the cluster centre, possessing a more extended distribution. 
The k+a types, meanwhile, are half-way between the k and 
e(c) types. 

Finally, in Figure 13, we also plot the luminosity dis- 
tribution of field galaxies at z 0.1 to compare with our 
composite cluster sample. Our field sample contains an ob- 
vious lack of bright k+a/a+k galaxies in comparison to the 
cluster sample. We display the combined bright k+a/a+k 
cluster galaxy spectra in Figure 2. Whilst these combined 
spectra do possess enhanced He and H7, there is a concern 
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r, / Mpc log„ (E) 

Figure 12. Median EW [On]A3727A as a function of projected radius from the cluster centre (left) and local galaxy density (right). 
The error bars are la on the median values. Each bin contains 100 cluster galaxies. The upper and lower curves are the 75tli and 25th 
quartiles of the cluster EW distribution. The horizontal solid lines are the 75th and 25tli quartiles of the field sample. The dotted vertical 
line denotes the critical radius (density) obtained by Gomez et al. (2003). There is a break in cluster EW [Oii]A3727A at about 3-4 Mpc 
(~ few galaxies per Mpc^) which corresponds well with the Gomez et al. (2003) and Lewis et al. (2002) results. 




Figure 13. Luminosity distributions of the different spectral types in the combined cluster (solid line) and the field sample (dotted line). 
The cluster k types are much brighter than the others; indeed the relatively flat k+a/a+k distribution may be due to mixing bright 
k-types with real, fainter k+a/a+k's. Meanwhile, the field sample appears to have a deficit of bright k+a/a+k galaxies in compajrison to 
the cluster luminosity distributions. 
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f / Tgoo 



Figure 14. Velocity-radius distributions of the different spectral types in the combined cluster. For reference, we superimpose the mass 
model from Figure 5 and a line at Av/ct^ = 0. The more active types possess a larger velocity dispersion and are more likely to reside 
in the outer regions of the clusters. 



over whether the bright k+a/a+k population seen in the 
clusters axe 'real' E+A galaxies. For example, Blake et al. 
(2004) find that of their 243 H5 strong E+A galaxies, 60 per 
cent have indications of ongoing star formation (Ha emis- 
sion; H7 and H/3 arc presumably subject to emission-filling 
and [Oil] is dust-obscured). Our cluster members do not 
have sufficient spectral coverage to reach Ha, so this can- 
not be checked. But in the absence of significant [Oil], it 
could be that like their E-l-A sample, ours may also be dust 
affected (e.g. dusty disks). 

We cam assess what would happen to our sample if we 
had have used a non-MORPHS, stricter definition of E-t-A 
galaxies: EW's of H5, H7 and H/3 all being greater than 4 A. 
The result is a loss of over two thirds of our k-|-a/a-|-k sample 
and we end up with some rather strange 'k' types that have 
very strong HJ EW! After all, there is much scatter in the 
H(5 versus H7 EW relationship (see Figure 1 of Blake et al. 
2004). Moreover, if k+a galaxies are selected on this strict 
basis, it is very unlikely that they will have any Ha emission 
(Blake et al. 2004). 

To test the degree of misclassification in the sample, we 
now perturb the EW measurements of [OH] and H(5. This 
is done by adding a random Gaussian sampling of the EW 
errors to the original measurements 100 times. With each 
new measurement, we compute what spectral class would 
have been assigned to the galaxy. In the combined cluster 
(Table 4) this re-typing would increase the number of k type 
galaxies at the expense of the k-|-a types. However, this in- 
crease is not very significant: no more than la. In the com- 



bined field fractions (Table 4), the perturbations result in an 
increase of the fractions of k and k-|-a types at the expense 
of a few e(c) and e(a) types. Again, however, this change is 
only by ~ 1.5cr. 



5.6 Backsplash 

We have already noted that some of our results may be 
considered to be due to galaxy mixing from 'backsplash' 
galaxies. Backsplash galaxies are those that have already 
been inside the cluster virial radius ^200) but are now 
observed outside this radius (e.g. Gill, Knebe & Gibson 
2005; Rines et al. 2005; Mamon et al. 2004; Balogh et al. 
2000). Gill et al. (2005) show through their simulations that 
the backsplash population can be detected in a straight 
forward fashion as they will exhibit distinct kinematics: a 
more centrally peaked velocity distribution in the interval 
1.4 < rj,/r2oo < 2.8. A non-zero velocity distribution peak 
would indicate that the galaxies in this interval are likely 
on their first infall. Rines et al. (2005) use this interval to 
show that backsplash cannot be wholly (or even primarily) 
the driving force behind galaxy transformations in clusters. 

We are also able to test the backsplash scenario through 
our observations (Figure 14). A histogram of the galaxy ve- 
locities in the interval 1.4 < Tylr^m < 2.8 is presented in 
Figure 16. This velocity distribution is not peaked at zero 
or at a larger value - therefore our cluster population at 
these radii is neither a pure backsplash population (which 
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Figure 15. Cumulative radial fraction per galaxy type in the 
combined cluster from Figure 14. The median radii for each type 
is noted in Table 6. The k-types are more centrally concentrated 
than the other types with poststarforming k+a types being half- 
way between k and e(c) types. 

would be expected to peak at zero velocity) or purely in- 
falling for the first time (which would be expected to peak 
at about 0.3-0.5 \/S.v\/(Jz in Figure 16; approximately 400 
kms~^; Gill et al. 2005). This mixture of both backsplash 
and infalling galaxies is in line with the recent results of 
Rines et al. (2005). Further, our data show no significant 
difference in velocity distribution for emission line galajcies 
(c(a) and e(c)) versus non-emission line galaxies (k, k+a 
and a-|-k) using a standard K-S test. Therefore, we comple- 
ment the result of Rines et al. (2005) that some emission line 
galaxies are likely to be backsplash ones and a pure back- 
splash model cannot account for the data. This is inline with 
the results of Gomez et al. (2003) and Lewis et al. (2002) 
that demonstrate that star formation is primarily a func- 
tion galaxy density. Our earlier suggestion that local galaxy 
density is a more fundamental parameter than radius from 
a cluster centre is also reinforced because of this. 



6 CONCLUSIONS 

In conclusion, this work presents, and analyzes, the spectro- 
scopic observations of rich X-ray luminous galaxy clusters 
at 2 ~ 0.1. In total, we have amassed over 7500 moderate 
resolution spectra, giving an unprecedented view of galaxy 
clusters from their core out to beyond their virial radii. Of 
these spectra, about 1600 are galaxy cluster members across 
eleven z ~ 0.1 clusters with good signal to noise ratios. Eas- 
ily two thirds of these cluster members are outside r2oo, 
making our survey one of the widest-field, high complete- 
ness surveys of clusters at these redshifts to date. 

• We confirm the result of P02 (see also Wake et al. 




Figure 16. Distribution of cluster galaxies velocities in the inter- 
val 1.4 < rp/r200 < 2.8 as derived from Figure 14. Our galaxies 
are not significantly peaked around zero velocity, nor at a larger 
velocity, indicating that LARCS has a mixed population of back- 
splash galaxies and those that are infalling into the cluster for the 
first time. 



2005) of a radial colour gradient at fixed luminosity of 
d{B - R) /drp = -0.011 ± 0.003 from the cluster core out to 
at least 6 Mpc. The composite cluster also displays a colour 
gradient with local galaxy density of d{B — R) / dlogiY,) = 
—0.062 ± 0.009 across three orders of magnitude in local 
galaxy density. In addition to these photometric gradients 
detailed here and in P02, we have also found spectroscopic 
gradients across our clusters, indicating the luminosity- 
weighted mean stellar population age in galaxies decreases 
with increasing radius [dRS/drp — 0.35±0.06), and decreas- 
ing local galaxy density (dH(J/d/ofl(E) = 0.66 ± 0.15). The 
effect is slightly stronger with local galaxy density than with 
radius. 

• A break in the [Oii]A3727A EW versus radius relation 

occurs at roughly rp = 4 Mpc. This is equivalent to ~ 1-few 
galaxies per Mpc'^ and is consistent with the result found by 
others (Lewis et al. 2002; Gomez et al. 2003). Such a large 
radius is somewhat surprising for processes which rely on 
high gas or galaxy densities to influence the starformation 
of cluster galaxies. We therefore suggest that effects on the 
star formation histories of galaxies take place in the infalling 
substructures, wliicli can have densities in the star forma- 
tion break range of ~few galaxies per Mpc^, in a manner 
independent of the cluster environment. This break in ra- 
dius is significantly different to the study made by Moran 
et al. (2005) who find a break at a radius of ~ 1 Mpc. This 
may be due to selection effects (Moran et al. 2005 sample 
galaxies down to M* + 1.6, have a smaller radial extent and 
define 'field' values differently) or that they examine only a 
single cluster (the cluster is at z = 0.4 and has an unusually 
large scatter in fundamental plane residuals in the cluster 
core region). 

• Even though our clusters are selected above a fixed 
X-ray luminosity, the cluster themselves show significant 
differences between their galaxy populations (galaxies with 
old stellar populations; poststarburst galaxies and currently 
starforming galaxies). The ratio of actively starforming to 
passive types in our clusters spans a large range 1 : 1.5 
through 1 : 9.5. There appear to be no clear correlations 
between the fractions of different spectral types and global 
cluster properties such as Lx, (Jz or (coarse) cluster mor- 
phology. However, we do find a mild anti-correlation between 
the Gaussianity of clusters velocity distributions and their 
fraction of k type galaxies. We suggest that the differences 
between the populations are probably the product of their 
recent accretion histories and are driven by local processes. 
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• There exist populations of galaxies both blueward and 
redward of the CMR that have similar spectral characteris- 
tics (i.e. old stellar populations with little or mild on-going 
starformation) to those galaxies lying on the CMR. We spec- 
ulate that the redder population, which arc morphologically 
indistinguishable from the CMR population, are dusty k 
types. The bluer population, meanwhile, are significantly 
less centrally concentrated and therefore we suggest that 
they may be the low redshift analogues to the passively 
evolving spirals found at higher rcdshifts (e.g. Balogh ct 
al. 2002; Poggianti et al. 1999). These bluer galaxies may 
ultimately evolve on to the CMR. 

• By constraining our 2: ~ 0.1 sample to the same radial 
extent and absolute magnitude range as the MORPHs sam- 
ple, we find that the luminosity distributions and spectral 
populations of cluster galaxies have evolved since z ~ 0.4. In 
particular, our z ~ 0.1 sample has significantly (> 3cr) more 
k type galaxies than their higher redshift counterparts. We 
suggest that the highly starforming e(b) class has declined 
(although the significance of this difference is only marginal 
at ^ 2.5cr). Certain other spectral types (k-|-a, a-|-k, e(a) and 
c(c)) show no significant changes. 

• Our data do not highlight any single mechanism or sce- 
nario that can wholly account for the results presented here. 
For example, our population at 1.4 < rp/r2oo < 2.8 consists 
of a mixture of both backsplash galaxies and galaxies that 
are infalling for the first time. Therefore star formation in 
these systems must be more fundamentally a function of lo- 
cal galaxy density rather than radius from the cluster centre. 
The break in EW measurements would not favour mecha- 
nisms that rely on high gas or galaxy density, so a mech- 
anism such as strangulation cannot wholly account for our 
data. Mechanisms that only affect the high density cluster 
cores are not favoured either (e.g. ram pressure stripping or 
galaxy-galaxy harrassment). 

This work follows Pimbblet et al. (2001; 2002) and is 
the third paper in a series based upon the LARCS survey. 
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APPENDIX A: SPECTROSCOPIC 
CATALOGUES 

Tabic 7 gives the LARCS speetroseopie cluster catalogues 
for the observations used in this work. The full version of 
our spectoscopic catalogues will be made available in Syn- 
ergy, the on-line version of the Monthy Notices of the Royal 
Astronomical Society. We note that a sub-set of these data 
have already been used by Pimbblet et al. (2005) and Krick, 
Bernstein & Pimbblet (2006). 



APPENDIX B: OTHER STRUCTURES 

As can be seen from Figure 4, there arc a rmrnbcr of other 
ovordense structures (i.e. walls, clusters, etc.) along the line 
of sight to our primary LARCS cluster targets. Here, we at- 
tempt to characterize these structures by applying the ZHG 
technique to them. In doing so, wo have only chosen those 
peaks in Figure 4 that arc significant (N(gal)>~20 in a red- 
shift range of dcz — 1000 kms^^) and relatively easy to 
isolate (i.e. no large spread in redshift space). As with the 
LARCS clusters, we make an initial guess of the velocity 
clipping to use by manual inspection of the velocity his- 
tograms (Figure 4). The results of this analysis are presented 
in Table 8. Note that since many of these structures are not 
wholly contained within the LARCS 2 degree field, we make 
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Table 7. LARCS spectroscopic cluster catalogues. Entries of NA indicate that whilst the R band was observed, the B band was not 
(Pimbblet 2001; POl). 



LARCS 


RA 


Dec. 


R Mag 


(B-R) 


cz 


Rtdr 


Identifier 


(J2000) 


(J2000) 






(km-l) 




a22_01:1148 


00 23 13.87 


-25 45 34.0 


17.40 ± 0.01 


0.39 ± 0.01 


-31 ± 16 


15.4 


a22_01:392 


00 24 51.43 


-25 50 57.2 


18.38 ± 0.02 


2.22 ± 0.07 


101831 ± 72 


2.6 


a22_02:1155 


00 21 32.76 


-25 44 10.9 


19.31 ± 0.05 


2.24 ± 0.18 


87177 ± 127 


7.2 


a22_02:1204 


00 22 9.26 


-25 45 19.8 


18.44 ± 0.02 


1.08 ± 0.03 


28267 ± 51 


5.6 


a22_02:1225 


00 21 59.53 


-25 45 29.9 


18.05 ± 0.02 


0.85 ± 0.03 


36457 ± 36 


4.8 


a22_02:1226 


00 21 54.37 


-25 45 29.9 


19.83 ± 0.05 


1.58 ± 0.08 


94392 ± 17 


0.1 


a22_02:125 


00 21 53.08 


-25 33 3.3 


17.22 ± 0.01 


0.65 ± 0.01 


67000 ± 37 


2.6 


a22_02:1258 


00 21 46.66 


-25 45 56.8 


19.74 ± 0.05 


1.62 ± 0.10 


34119 ± 19 


0.8 



Table 8. Other structures identified along the line of sight to the primary LARCS cluster targets. Entries of 'NED02' refer to objects 
that have multiple entries with the same name in the NED database - the 02 denotes that this is a secondary object not associated with 
the primary object of the same name. 



Structure 


cz Range 


cz 




N(gal) 


Conirnents 


Idcnl iii(*r 


(knis J-) 


(knis \) 


(knis J-) 






LARCS-22-Wall 


30000-35000 


33811±88 


6191^4 


>49 


Part of candidate superclustcr? (Pimbblet et al. 2005) 


Abell 1079 


38000-42000 


39671±64 


719+50 


46 


Attempt to localize Abell 1079 in field of Abell 1084 


LARCS- 1437- A 


24500-27500 


25874±90 




40 


Possibly Abell 1459(?) 


LARCS-1650-A 


33000-38000 


35981±99 


"700+82 


54 




LARCS-1651-A 


40000-45000 


41541±223 


,.„„+193 


37 


Possibly Abell 1658(?) 


LARCS-1651-B 


53000-57000 


54972±116 




48 


ZwCl 1300.9-0244(?) 


LARCS-2055-A 


8500-12000 


10366±74 


495tg 


44 


Abell 2052 


LARCS-2055-B 


22000-26000 


22766 ± 152 


sei^if 


31 




LARCS-2055-C 


51000-57000 


53825 ± 174 


055+155 


29 




LARCS-2104-A 


16000-20000 


17567 ± 184 




23 


Abell 2103(?) 


LARCS-2204-A 


16000-18500 


17441 ± 50 


3601^2 


51 




LARCS-2204-B 


22000-25000 


23810 ± 41 


3931^^ 


91 




LARCS-2204-C 


86000-88000 


87721 ± 36 


235 ^0 


38 




LARCS-3888-A 


21000-24000 


22691 ± 118 


781+10" 


43 


Possibly Abell 3896 NED02(?) (Batuski et al. 1999) 


LARCS-3888-B 


58000-65000 


61232 ± 224 


1449+1«^ 


41 




LARCS-3921-A 


37500-42200 


40208 ± 112 


945+?i 


70 


Abell 3921 NED02 (Katgert et al. 1996) 


LARCS-3921-B 


47000-49000 


48114 ± 99 


476l«^ 


22 




LARCS-3921-C 


60000-65000 


62893 ± 145 


948+gf 


42 


AM 2250-633 


LARCS-3921-D 


85000-90000 


87505 ± 170 


1050ll*oI 


37 





no attempt to find the centroid of them. However, we do 
use individual members positions to aid in identifying hkely 
corresponding clusters from the Uterature. 
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